Practical and accurate calculations of Askaryan radiation. 
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An in-depth characterization of coherent radio Cherenkov pulses from particle showers in dense 
dielectric media, referred to as the Askaryan effect, is presented. The time-domain calculation 
developed in this article is based on a form factor to account for the lateral dimensions of the 
shower. It is computationally efficient and able to reproduce the results of detailed particle shower 
simulations with high fidelity in most regions of practical interest including Fresnel effects due to the 
longitudinal development of the shower. In addition, an intuitive interpretation of the characteristics 
of the Askaryan pulse is provided. We expect our approach to benefit the analysis of radio pulses 
in experiments exploiting the radio technique. 
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I. INTRODUCTION 

In 1962 Askaryan [l| proposed to detect Ultra-High 
Energy (UHE) cosmic rays and neutrinos observing 
the coherent radio pulse from the excess of electrons 
in a shower developing in a dense dielectric and non- 
absorptive medium. The scaling of the emitted power 
with the square of the particle energy, which has been 
experimentally confirmed in accelerators , makes the 
technique very promising for the detection of UHE par- 
ticles and has motivated a variety of past and present 
experiments along with some in the planning 

stages 

Key to the success of these initiatives is an accurate 
and computationally efficient calculation of the radio 
emission properties due to the Askaryan effect in UHE 
showers. The problem of computing the coherent radia- 
tion from particle showers can be approached in a variety 
of ways. Purely Monte Carlo methods have been devel- 
oped to simulate the induced showers in dense media. 
One can obtain the contribution to the radiation from 
every particle track in the shower from first principles 
(Maxwell's equations), and add the contributions, which 
automatically takes coherent effects into account. This 
approach has been applied to calculate the Fourier com- 
ponents of the radiation (i.e. in the frequency-domain) 
|16l - l28| . and only recently to the calculation of the ra- 
diation as a function of time (i.e. in the time-domain) 
[29j . Similar methods have also been applied to the cal- 
culation of radio emission in atmospheric showers (sol - fssj 
in which the Askaryan effect is not the dominant mech- 
anism and will not be addressed in this paper. Monte 
Carlo methods have the advantage that the full complex- 



ity of shower phenomena is accounted for, the influence 
of shower-to-shower fluctuations can be addressed, and 
the dependence on the type of primary particle, hadronic 
model, along with any other assumptions can be stud- 
ied with high accuracy. However, purely Monte Carlo 
methods are typically very time-consuming, especially at 
ultra- high energies and approximations are required [27| . 

Another numerical approach currently being developed 
is the application of finite difference in the time-domain 
(FDTD) techniques [s^. The idea is to discretize space- 
time and propagate the electric and magnetic fields by 
approximation of Maxwell's differential equations into 
difference equations [s^l ■ FDTD techniques have the ad- 
vantage that they can be easily adapted to computing 
the effects of dielectric boundaries and index of refraction 
gradients and can be linked to an accurate Monte Carlo 
simulation of showers in dense media. The FDTD tech- 
nique is however rather computationally intensive (36j . 

Analytical approaches have also been developed. In 
these methods the charge development in the shower is 
approximated as a current density vector [s^ . Typically, 
parameterizations of the longitudinal and lateral profile 
of the showers are used to describe the main features of 
the space-time evolution of the charge distribution. In 
this approach one calculates the vector potential by inte- 
grating the Green's function to obtain the electric field. 
These integrals are in general difficult both numerically 
and analytically, but they can be greatly simplified by the 
use of approximations [38|. These methods are usually 
less time-consuming but also less accurate than purely 
Monte Carlo simulations. One limitation is that show- 
ers elongated due to the Landau-Pomeranchuk-Migdal 
(LPM) effect [H, cannot be parameterized easily due 
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to large shower-to-shower fluctuations and the emitted 
radiation is known to depend stron gly on the particu- 
lar longitudinal profile of the shower |18l-l20l[. Analytical 
techniques have also been applied with different levels 
of sophistication to the calculation of radio emission in 
atmospheric showers j4ll - l43 |. 

It is clearly important to be able to calculate coher- 
ent radiation due to the Askaryan as well as other ef- 
fects with a variety of techniques as each has its own 
set of advantages and disadvantages. Our goal in this 
work is to provide a calculation method that is both fast 
and able to reproduce all the essential characteristics re- 
sulting from detailed Monte Carlo shower simulations. 
Semi-analytical techniques are a very good option. The 
general idea behind this method is to obtain the charge 
distribution from detailed Monte Carlo simulations as 
the input for an analytical calculation of the radio pulse. 
Other semi- analytical methods in dense media [T8l - [20l [29j 
and in the atmosphere 45,] have been attempted in the 
past. In dense media, the frequency spectrum of the ra- 
dio emission due to the Askaryan effect has been shown 
to be easily obtained from the Fourier-transform of the 
longitudinal profile of the shower [20j . This technique re- 
produced the frequency spectrum as predicted in Monte 
Carlo simulations with a high degree of accuracy, but 
only for angles away from the Cherenkov angle in the 
Fraunhofer approximation. Complementary, in the time- 
domain, it was shown that the electric field away from the 
Cherenkov angle and in the far-field regime can be very 
accurately calculated from the time-derivative of the sim- 
ulated longitudinal development of the excess charge [2^ 
(see also |42|). 

This paper presents a semi-analytical calculation that 
is able to reproduce the electric field in the time-domain 
at all angles with respect to the shower axis in both the 
far-field (Fraunhofer) and "near-field" Q regions of the 
shower when compared to a full Monte Carlo simulation 
such as the well-known Zas-Halzen-Stanev (ZHS) code 
[TtI [29l . |46| . The technique is computationally efficient 
since it only requires the convolution of the longitudinal 
charge excess profile with a parameterized form factor to 
fully reproduce the coherent radiation effects from par- 
ticle showers in a homogeneous dielectric medium. Once 
the longitudinal shower profile is obtained, the electric 
field can be calculated with a simple numerical integral. 
It is worth remarking here that the longitudinal devel- 
opment of extremely energetic showers can be obtained 
quickly and with high precision, using hybrid simulation 
techniques which consist on following only the highest 
energy particles in the shower while accounting for the 
lowest energy particles with parameterizations. In partic- 
ular, the complexity of the longitudinal profile of showers 



By "near-field" we mean a region in which the Fraunhofer ap- 
proximation is not vaUd because of the longitudinal dimensions 
of the shower and not the region where the Coulomb field asso- 
ciated with the charge excess cannot be considered negligible. 



affected by the LPM effect can be very well reproduced 
withjrybrid techniques [l8l - l20l | (for an example see Fig. 5 
in 

The semi-analytical method described in this work is 
well suited to obtain the time-domain radio emission 
due to electromagnetic showers, for all observation an- 
gles both in the far-field and the near-field regions of 
the shower. The approach can be used in practically all 
experimental situations of interest since it only begins to 
show significant discrepancies when the observer is at dis- 
tances comparable to the lateral dimensions of the shower 
(< 1 m in ice). Since the typical distance between an- 
tennas in experiments such as the Askaryan Radio Array 
(ARA) 14] is ^ 10 — 100 m, we expect the results to be 
accurate enough in most practical situations. 

We expect our results to benefit experiments exploiting 
the radio technique. They can be used in detector simu- 
lations to test the efficiency for pulses observed from var- 
ious directions. In particular, with our approach one can 
test the ability to detect the craggy pulses resulting from 
the LPM effect. The calculation can also be implemented 
in the data analysis of experiments by using likelihood 
functions aimed at the reconstruction of the longitudinal 
charge excess profile from a detected pulse. If the longi- 
tudinal distribution is consistent with the elongation and 
multiple-peaked structure due to the LPM effect, this can 
be used for neutrino flavor identification since it is only 
expected in UHE showers due to electron neutrinos. 



II. MODELING ASKARYAN RADIATION 

The case we are interested in this paper is the radiation 
due to the charge excess of a shower in a linear dielectric 
medium such as ice, salt, or silica sand. We use SI units 
all throughout this work. The Green's function solutions 
to Maxwell's equations provide the potentials $ and A 
given a charge distribution p with current density vector 
J = pv. Assuming a dielectric constant e and magnetic 
constant p the solutions in the Coulomb gauge (V-A = 0) 
can be written as 

rav (1) 

47re Jx-x'l 

A(x,t) = £ £i(^^(^|x_x'l-(i-0)rfVdt' 

(2) 

where the delta function gives the observer's time t de- 
layed with respect to the source time t' by the time it 
takes light to reach the observation point x from the 
source position at x'. The transverse current is given 
by = — u X (u X J) where u = (x — x')/]x — x'] is 
the unit vector pointing from the source to the observer. 
The non-trivial proof that the Jj^ above is the only com- 
ponent relevant to the radiation part of the field is given 
in 113. 
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Radiation calculations are typically performed in the 
Lorentz gauge (V-A + n^ c-^d<^/dt = 0) [M HI El Ell , 
with n the refractive index of the medium. However, our 
primary interest is to derive an approach that can be 
easily applied to numerical radiation calculations. In the 
Coulomb gauge, the scalar potential only describes near- 
field terms which can be ignored for our purposes. This 
simplifies the computation of the radiative electric field 
E = — V$ — dA/dt to a simple time derivative E = 
—dA/dt. Thus, all that is needed is a computation of 
the vector potential as a function of time at the position 
of interest. 

The radiation of a particle shower in a dense medium 
is obtained by treating it as a current density J with its 
main features depicted in Fig[TJ The evolution in space- 
time of the excess charge in a shower can be modeled 
as a pancake 6(z' — vt') traveling with velocity v along 
the z-axis. The net charge profile of the shower Q{z') 
rises and falls along the shower direction z' and spreads 
laterally in x' and y' . The velocity vector v is primarily 
directed in the shower axis direction z but may have a 
small lateral component and a lateral dependence due to 
scattering of particles in the shower. This may seem like 
an unnecessary complication but the small scatter will 
lead to an observable asymmetry of the Askaryan pulse 
in the time-domain. The speed v is assumed to be close 
to the speed of light for particle showers of interest. The 
associated current density vector can be modeled as a 
cylindrically symmetric function 

J(x', t') = v(r', 0', z')f{r', z')Q{z')diz' - vt') (3) 



where r' = x'^ + y'^ is the cylindrical radius. The func- 
tion /(r', z') represents the lateral charge distribution in 
a plane transverse to z' depicted in the bottom of Fig. [1] 
This model of the current density is similar to the ones 
used in [2^ except that v is allowed to have first order 
radial components. 

The vector potential is then given by 



A(x,t) = /°° dt' r dz'Q{z')5{z' ~ vt') 

J —CO J —oo 



oo p27r 



dr'r' I d4)'f{r',z')v^{r',(t)',z') 
10 Jo 
6{n\x-K'\/c^{t-t')) 



(4) 



where 0' is the azimuthal angle in cylindrical coordinates 
and v_L(r', </>', z') generally depends on r' , (f>' , and z' . 

In its full glory Eq. (|4|) seems rather intractable. The 
lateral distribution f{r',z') is the most difficult part to 
model. It is for this reason that the semi-analytical mod- 
els developed in [23| and [29^ ignored the lateral devel- 
opment of the shower, and hence were not able to accu- 
rately describe Askaryan radiation near the Cherenkov 
angle where the lateral distribution is known to deter- 
mine the degree of coherence of the emitted radiation 
[l7| . Modeling the radiation due to the lateral distribu- 
tion of the shower can be attempted using the standard 
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FIG. 1. Geometry of a high energy particle shower. The 
top figure shows a side view {r — z plane) of the charge excess 
profile. The shower front propagates at velocity v and is mod- 
eled as a thin pancake S{z' — vt'). The charge evolution of the 
shower front traces out the longitudinal profile Q{z') shown 
here as an asymmetric distribution. The lateral spread of the 
shower due to Coulomb scattering, shown in (blue) arrows, 
is modeled by a velocity vector v — v(r',(^') that is mostly 
directed along the shower axis with small radial components. 
The scatter results in observation angles that differ from the 
nominal angle 9^ relative to the shower axis and will lead to an 
asymmetric pulse. The radiative portion of the electric field 
lies along the vector p which is the orthogonal projection of 
z along the direction u from the source to the observer. The 
bottom panel shows a frontal view (r — <^ plane) of the shower. 
The lateral charge excess distribution f{r') is represented in 
(blue) shading. 



NKG function |49|]. However, this parameterization has 
singularities that make the integrals particularly difficult 
to solve and interpret, and it assumes constant particle 
velocities parallel to the shower axis which is not suffi- 
ciently accurate as will be shown further below. 

In the following we will show that the radiation due 
to the lateral distribution of the shower can be parame- 
terized, and this parameterization can be used to predict 
the emitted radiation matching that obtained from the 
full ZHS simulation [17, 29] with great accuracy. More- 
over, the lateral distribution of the shower is mainly due 
to low energy processes such as Coulomb scattering, and 
as result the shape of the parameterized radiation is inde- 
pendent of shower energy in the energy range of interest 
as will be shown below. On the other hand, the longitudi- 
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nal distribution Q{z') can change dramatically depending 
on the energy of the shower due to the LPM effect, and 
needs to be obtained in a Monte Carlo simulations such as 
ZHS. We will also show that our calculation of Askaryan 
radiation works in both the near and far-field approxi- 
mations, as long as the lateral coordinates are treated in 
the far-field, i.e. the observation of the shower occurs at 
a distance larger than its lateral dimensions {'^ 1 m in 
ice). 



A. The vector potential at the Cherenkov angle 

Let us first consider the Fraunhofer approximation for 
radiation emitted by a shower. This implies expanding 
|x - x'l in Eq. ^ as 



X 



(5) 



where R = |x|. For an observer looking in the direction 
u — x/|x| = (sin 6* cos 0, sin 6* sin 0, cos 0) in spherical co- 
ordinates, and assuming without loss of generality that 
(/) = 0, the above expansion can be written as 



R 



r sm tf cos i 



(6) 



Approximating the denominator of the vector potential 
in Eq. @ by |x — x'| « R but keeping the approximation 
in Eq. ([6]) in the argument of the 5- function we get. 



A{0,t) = 



OO /'OO 

dt' 



dz'Q{z')5{z' - vt') 



AttR 

dr'r' I d(j>'f{r',z')v^{r',4)',z') 



(7) 
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z' cos ( 



Integrating over the source time t' results in 
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where v_L(r', </>', z') = — u x u x v(r', </>', z') is the trans- 
verse projection of the unit velocity vector. 

In our model we make the assumption that the shape of 
the lateral density and the particle velocity depend only 
very weakly on z', i.e. f{r', z') ~ f{r') and v(r', 0', z') ~ 
v(r', (/)'). With these approximations Eq. ([8|) can be writ- 
ten as: 
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(9) 
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FIG. 2. The vector potential from electromagnetic showers 
in homogeneous ice (density p — 0.924 g cm~^ and refractive 
index n — 1.78, 9c ~ 55.8°) observed at the Cherenkov angle 
for various energies from the ZHS simulation. The functional 
behavior is identical except for an overall scaling factor that 
is directly proportional to the shower energy. The width is 
determined by the lateral distribution of the shower while 
the asymmetry is mainly due to the radial spread of particle 
tracks; both are the result of the Coulomb scattering in the 
medium. 



where we have defined the function F as: 

F (^t 
1 

sin 9 



nR 1 


"1 


c 


V 


/"OO 




I "••'••'J 


1 








,^vx(r',0') 



(10) 



nr sm t) cos ( 



nR 
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Note that we have explicitly included a factor sin 6* in 
Eq. ([9]) for convenience, anticipating that the radiation 
is mainly polarized in the direction transverse to the ob- 
server's direction. 

The vector function F contains the radial r' and az- 
imuthal cfi' integrals and can be considered as an effec- 
tive form factor that accounts for the lateral distribution 
of the charged current density, quite analogous to that 
obtained in the frequency domain in (ssj . 

At the Cherenkov angle we have 1/v — n cos 9c /c = 
and Eq. ([9|) results in: 



dz'Q{z') 

(11) 
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where F at the Cherenkov angle is obtained from 
Eq. (UnD: 



Fit-— U-V / dr'r'fir 

c J smfc'c Jo Jo 

I nR nr' sin Oc cos </>' 

I t 
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(12) 



Using symmetry arguments we can now project the 
form factor F in only two orthogonal directions using 
unit vectors u along the observation direction and in the 
direction of p = — u x (u x z): 



F = Fp p + i^„u 



(13) 



The direction of p has been chosen perpendicular to the 
observation direction and lying on the plane defined by z 
and u as shown in Fig. [J]. This has been done anticipat- 
ing the expected polarization of the radiation mainly in 
the direction of p in order to make the orthogonal com- 
ponent along the direction of observation negligible. 

The form factor F is a medium dependent function 
that accounts in an effective way for the radial and az- 
imuthal interference effects due to the lateral structure 
of the shower including possible directional variations in 
the velocity vector v of the shower particles. In principle 
F can be obtained from analytical solutions of the cas- 
cade equations, but this is rather involved. Alternatively 
one could use standard parameterizations of the lateral 
distribution function of high energy showers such as the 
NKG d^. However, as stated before, this does not ac- 
count for the radial components of the velocity and gives 
results symmetric in time which are qualitatively differ- 
ent from results obtained in detailed simulations such as 
that shown in Fig. [5] The trick is to extract F from 
simulations that effectively account for the radial com- 
ponent of particle velocities and the lateral distribution 
of the excess charge which are responsible for the time 
asymmetry characteristic of simulations. The basic idea 
behind this article is that the form factor is obtained 
from the vector potential at the Cherenkov angle and, as 
we will show later, the emission at other angles is easily 
related to that at the Cherenkov angle. 

The vector potential at the Cherenkov angle in the 
time domain is calculated with a detailed shower simula- 
tion and parameterized for practical purposes. The form 
factor F can be obtained directly equating Eq. (fTTj) to the 
vector potential as obtained in the simulation. As antici- 
pated the Fu component is typically below 1 % of Fp and 
it can be neglected, and Fp can be simply obtained from 
the following equation: 



AttR 



sm Uc 



F„ t 



nR 



LQtot P (14) 



^ By use of vector identities the direction of p can be shown to be 
that of z — (z ■ u)u which is more apparent. 



where LQtot = J dz'Q{z') has been referred to as the 
excess projected track- length [l^. Taking the absolute 
value of Eq. ([T4l) the functional form of Fp is given by: 



F„\t- 



nR 
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sm Oc 



(15) 



where A{9c,t) = |A(0c',OI- Pp represents the average 
vector potential at the Cherenkov angle per unit excess 
track length - given by LQtot - scaled with the factor 
At:R/pl. 

Detailed simulations of electromagnetic showers per- 
formed with the ZHS Monte Carlo code in ice produce a 
consistent time-dependent vector potential at all energies 
of interest as shown in Fig. [2] The results for homoge- 
neous ice can be parameterized by: 



RA{ec,t) = -4.5 X 10"^* [V s 



E 



[TcV] 



exp 



exp 



(-a^)+(l + 2-87|tr3 ifi>0 (16) 



0.030 



(l + S.OSltl)-^-^ ifi<0 



where E is the energy of the shower in TeV and t is the 
observer time in ns. The result is accurate to within 5%. 
The shape of A{9cjt) depends very weakly on shower 
energy, while the normalization is proportional to the 
energy as becomes evident in Fig. [21 Note also that at 
higher energies the fluctuations are reduced because the 
number of particle tracks increases almost linearly with 
shower energy. 



B. The radiation in the far-field 

Given the time domain parametrization of the radi- 
ation at the Cherenkov angle in Eq. (ITB)) we will first 
obtain the pulse as seen by an observer in the far-field 
at any observation angle. The integral in Eq. ^ can be 
written as 
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dz'Q{z') 



F„ 



nR , 

t z' 
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n cos 6 



(17) 



This is a straightforward numerical integration with Fp 
given by Eqs. (|15|) and and the longitudinal distribu- 
tion of the excess charge Q{z') which can be obtained in 
a shower simulation or from a standard parameterization 
of the depth development if the LPM effect is absent. 
The corresponding radiative electric field is obtained by 
simply taking the derivative of the vector potential with 
respect to time E = —dA/dt. 

Eq. (IT71) tells us that the vector potential in the far- 
field region of the shower can be obtained as a convo- 
lution of the form factor Fp, which parameterizes the 
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FIG. 3. Results of the ZHS Monte Carlo compared to our cal- 
culation Eq. p7p for the radiation due to a 3 x 10^* eV elec- 
tromagnetic shower observed at 9 = 6c ~ 0.3° in the far-field 
in ice. The method is applied to the charge excess longitu- 
dinal profile obtained from the same simulation. Top panel: 
Vector potential in the time-domain with 10 ps time sam- 
pling (corresponding to a sampling frequency of 100 GHz) as 
obtained in ZHS simulations. This is compared to the calcula- 
tion presented in this work. The longitudinal charge profile is 
shown as the ID model presented in [29|] where the depth and 
charge are linearly rescaled to give the observer time and vec- 
tor potential. Middle panel: Electric field in the time-domain 
comparing our method to ZHS results. Bottom panel: The 
electric field amplitude spectrum obtained from the ZHS sim- 
ulation and the Fourier transform amplitudes of the electric 
field obtained from the calculation presented in this work. 
Note that the discrepancy between the time-domain electric 
field of the ZHS simulations and our results are due to the 
incoherent radiation at high frequencies. 



emission from the lateral distribution of the shower, with 
the longitudinal profile of the excess charge. The form 
factor Fp is a function that has to be evaluated at the 
time t at which the observer in the far field sees the por- 
tion of the shower corresponding to the depth z'. That 
time is given hy t = nR/c + z'/v — z'ncosO / c. We 
have made the only assumption that the shape of Fp de- 
pends weakly on the stage of longitudinal evolution of the 
shower. At the Cherenkov angle, the far-field observer 
sees the whole longitudinal development of the shower 
at once i.e. z'/v = z'ncosBc /c in which case Eq. (fT7|) 
reduces to the vector potential at the Cherenkov angle 
given by Eq. ([14]). 

In Fig. 13] we show an example of the vector potential 
and electric field in the time-domain due to an electro- 



magnetic shower with energy E = IQQ EeV from the ZHS 
simulation. The fields are observed in the Fraunhofer re- 
gion at an angle 6 = — 0.3° and they are compared to 
our results obtained with Eq. P7)) using the longitudinal 
distribution Q{z') from the same simulation. The agree- 
ment between the vector potential obtained directly in 
the Monte Carlo simulation and the prediction of Eq. pT)) 
lies within a few percent difference in the region relevant 
to the pulse (top panel of Fig. [3]) . The difference between 
this calculation and the ZHS electric field in the time do- 
main is greater (middle panel of Fig. [3]), but as shown in 
the bottom of Fig. [3] this is due mostly to the incoherent 
emission of the shower at high frequencies. The Fourier- 
transformed amplitudes of the time-domain electric field 
obtained in our approach are based on smooth param- 
eterizations, while the frequency spectrum obtained di- 
rectly in the ZHS simulation includes incoherence effects 
coming from the fine structure of the shower at the indi- 
vidual particle level. 

Note that when f{r',z') — S{r')/r', i.e. if we neglect 
the lateral distribution of the shower, Eq. PT)) reduces to 
the 1-dimensional model in [2^ which fails at describing 
the features of the radio emission for angles close to the 
Cherenkov angle. To illustrate this, we show in the top 
panel of Fig. [3]the vector potential obtained in the one di- 
mensional (ID) model in j29j . which is a linear rescaling 
of the longitudinal charge excess profile, showing a clear 
disagreement with the results of the full ZHS simulation 
as expected. 



C. Askaryan pulses in the "near-field" 

We can now generalize Eq. (I17|) for an observer in the 
"near-field" region of the shower. In this case it is more 
natural to work in cylindrical coordinates and place the 
observer at (r cos (j), f sin 0, z). Without loss of generality 
we can again assume the observer is at = giving 



'I = v^r2 + r'2 - 2rr' cos (j)' + {z - z'Y 



(18) 



In dense media, the lateral distribution is in the scale of 
centimeters, which means that for all practical purposes 
the observer is at any given instant in the far-field region 
with respect to the lateral distribution. The idea is to 
solve the vector potential in Eq. (JJ]) using the Fraunhofer 
approximation to account for the lateral distribution at 
any given time t' . We expand Eq. (|T5)) to first order in 
r' giving 



|x- x' I = v/?^^ + (r^^O^- r' sin 6l(z') cos 0'. (19) 



where sin0(z') = r/^Jr'^ -\- [z — z')-^, but we take into ac- 
count that the distance in the denominator of the vector 
potential depends on the time t' or equivalently on the 
position z' in the shower as \/r'^ + {z — z'Y . This is in 
contrast to the case of the far field calculation in which 
the distance in the denominator of the vector potential 
is constant and equal to R. 
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FIG. 4. Sketch illustrating the polarization dependence with 
time. As the shower evolves the observer sees the radiation 
contributions at different angles corresponding to different 
source positions z . This means that the observer will see 
the polarization vector changing as a function of time t. 



We proceed as in the case of the far-field calculation as- 
suming that f{r',z') ~ f{r') and v{r' , cj)' , z') ~ v(r', (/>'). 
After integrating over t' the vector potential in Eq. 1^ 
can be written as, 



A{r,z,t)^^ r dz'- 



yjr"^ + (z - z'Y 



-oo 



V 



z' n^J -f (z — z') 



^ — nn nr 



nr smf/ z cosi 



(20) 

where the transverse projection of the velocity vector 
v_L(r', 0', z') = — u(z') X [u{z') X v(r', (/)')] now intro- 
duces a new dependence on the longitudinal source co- 
ordinate z' due to the fact that in the near-field u(z') = 
{r+{z~z')z)/y/r'^ + {z — z')"^ depends on z'. If we define 
the form factor containing the radial r' and azimuthal (j)' 
integrals as in subsection II. A we obtain a similar expres- 
sion: 



^ z' n^r^ + (z — z'Y' 



III j^/ ti i\ ,0 ,z ) 



dr r \ d(j) f{r ) 
, Jo 



V 



g ( z' n^r^ + (z — z')2 ^ nr' sin 6'(z') cos 
\ V c c 



(21) 



Note that the form factor defined in Eq. (I^TI) has the 
same functional form as that defined in Eq. ([15]) and they 
only differ in the argument of the delta function. The 
form factor F as obtained in the far field can be applied to 
the near field (in relation to the longitudinal development 
of the shower) simply modifying its argument. 



Neglecting again the component of F parallel to u the 
vector potential in the near-field can be written as: 



A(r,z,t) 



47r 



Q{z') 

dz , p(z ) 



z' n^r^ + (z — z')^ 



(22) 



Note that a new z' dependence is introduced through 
the polarization vector p(z') = — u(z') x (u(z') x z). In 
Fig. 0] an example is shown of the z' dependence of p(z'). 
Also note that the sinf? term in Eq. pT|) has been ab- 
sorbed in Eq. ((22|) through p — sin Q p 

In the near field the explicit z' dependence is necessary 
because the longitudinal profile of the shower is observed 
at different angles. This means that parts of the shower 
observed at different depths will have differing polariza- 
tion vectors. This modification accounts exactly for the 
interference between different z' points along the shower 
development. This result has also been determined from 
a one dimensional current density model in (50| . 

Eq. fF2\ tells us that the vector potential in the near- 
field region of the shower can be obtained as a convolu- 
tion of the form factor - that parameterizes the in- 
terference effects due to the lateral distribution of the 
shower - and the longitudinal profile of the excess charge. 
The form factor function has to be evaluated at the 
time t at which an observer in the near field sees the por- 
tion of the shower corresponding to th e depth z^ Th at 
time is clearly given by t = z' jv ~\- n^r^ -I- (z — z')'^ /c. 
The difference between this expression and that obtained 
in the far field is that Fp is always evaluated at the time 
t an observer sees the position z' in the shower which is 
different for the far- and near-field regions. 

In Fig. [5] we show an example of the vector potential 
in the time domain for various observers in the near field 
region of the shower as obtained in full ZHS simulations. 
The shower has an energy E — 100 PeV and a longitu- 
dinal dimension of ^ 25 m, and the observer is placed 
at different positions z along the shower axis and at a 
fixed radial distance r = 10 m. The ZHS simulation also 
gives the longitudinal profile of the excess charge Q{z') 
which we have introduced into Eq. ([22]) to obtain the 
vector potential and compare it to that obtained directly 
by the Monte Carlo, also shown in Fig. [S] The agree- 
ment between the vector potential obtained directly in 
the Monte Carlo simulation and the calculation in our 
approach is remarkable. For distances to the shower axis 
larger than ~ 1 m the difference between our approach 
and the Monte Carlo is typically ~ 1 % or better for 
values down to ~ 3 orders of magnitude below the peak 
of the vector potential. This difference starts to increase 
gradually as the distance to the shower axis decreases 
and becomes comparable to the lateral dimensions of the 
shower, where the parameterization of the vector poten- 
tial at the Cherenkov angle given in Eq. (IT51) is not ex- 
pected to be valid. Since the typical distance between 
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FIG. 5. The Fresnel region vector potential for an electron 
induced 100 PeV shower in ice from the ZHS Monte Carlo 
compared to the results of Eq. H22|) . The longitudinal profile 
of the excess charge is convolved with the form factor Fp. The 
observers are located at {x, y, z) where z lies along the shower 
axis and y = Q. The vector potentials have been arbitrarily 
shifted in time by multiples of 50 ns for clarity. The spiked 
vector potential in the middle (blue) corresponds to an obser- 
vation where the peak of the longitudinal shower profile is near 
the Cherenkov angle. The vector potentials to the right of it 
(green and red) are due to observations above the Cherenkov 
angle while the vector potentials to the left (magenta and yel- 
low) are due to observations below the Cherenkov angle. Note 
the inverted order of the primary and secondary peaks. 



antennas in exp eriments such as the Askaryan Radio Ar- 
ray (ARA) [ij] is ~ 10 — 100 m, we expect our results to 
be accurate enough in most practical situations. 



III. THE APPARENT MOTION OF A CHARGE 
DISTRIBUTION. 

The temporal behavior of the vector potential traces 
the motion of a charged particle according to the retarded 



time. This is an old idea discussed by Feynman in [51 1 
applied to elucidate on the properties of synchrotron radi- 
ation. More recently, this approach has also been used in 
the Cherenkov radiation calculation due to linear tracks 
in the near field [s^ and one dimensional current densi- 
ties in [HO] • In this section we analyze the characteristics 
of our results in terms of the apparent motion of the 
charge density distribution to gain an intuitive under- 
standing of the radiation due to a particle shower devel- 
oping in a homogeneous dielectric medium. 
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FIG. 6. Illustration of the vector potential characteristics due 
to the apparent motion of the charge distribution for an ob- 
server seeing the region around shower maximum with angles 
close to the Cherenkov angle in ice (n — 1.78). In the bottom 
plot we show the longitudinal charge excess profile follow- 
ing a Greisen parameterization. The center plot shows the 
observer time as a function of the source time for a shower 
exceeding the speed of light in a medium. The solid (red) 
portion, corresponding to the minimum of the source and ob- 
server time relation, contributes significantly to the radiation 
near the Cherenkov peak. The resulting vector potential is 
shown in the right hand side. The Cherenkov radiation spike 
corresponds to the compressed mapping of the charge excess 
distribution to the observer time. 



A. Time delay effects 

The apparent motion of particles is encoded in the 
Green's function solution of the vector potential, Eq. (|4|) , 
by the argument of the delta function taking the source 
time t' to the observer time t: 



(23) 



The term |x — x'| traces the motion of the current density 
vector J at position x'(i') to determine the observer time. 

We can gain much insight into the properties of the 
vector potential resulting from the ZHS particle shower 
simulation, shown in Fig. [Sj by momentarily ignoring the 
lateral distribution of the shower. In this case the ob- 
server time t is given in terms of the source time t' by 



t^t' + 



(24) 



where we have substituted z' — vt' and r is the cylindrical 
radial position of the observer. 

When V < c/n, Eq. (j24p has a unique observer time 
corresponding to each source time. However, in the case 
of t; > c/n there always exist a range of observer positions 
such that Eq. has two source times corresponding to 
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every observer time (see an example in the top left panel 
of Fig. in]). In addition, a minimum value of t (different 
from the trivial minimum z'q corresponding the beginning 
of the shower at ip) may exist when n/3 > 1. In other 
words, the observer first sees the radiation corresponding 
to a depth in the shower z^j^ 7^ Zq ^^'^ then sees con- 
tributions from shower depths before and after arriving 
simultaneously. This minimum can be characterized by 
looking at the derivative of the retarded time relation, 



vt' 



+ {z- vt'Y 



(25) 



The extrema in the relation between the source time and 
the observer time are given by requiring the above equa- 
tion to be equal to zero. A solution exists only when 
V > c/n and indeed corresponds to a minimum value in 
the observer time imin- The corresponding shower co- 
ordinate z'- given by the source time t'- = z'-jv 

mm ^ ^ mm mm/ 

is 



(26) 



1 



The angle of observation Q with respect to the shower 
axis corresponding to the shower position z' is given by: 



cos Q 



{z-z') 



yjr"^ + {z- z' f 



(27) 



when z' — 



it is straightforward to show that the an- 



gle corresponds to the Cherenkov angle cos^c = l/('^/3), 
i.e. the minimum time imin at which the observer first 
sees the shower corresponds to the shower coordinate ly- 
ing at the Cherenkov angle. Its value is given by 



1 



(28) 



A peculiar consequence of these relations is that if an 
observer is placed at a position such that the shower is 
always seen with 6 < 6c then z'^-^^ corresponds to the end 
of the shower, which is an apparent violation of causality. 
In the case where > 9c then z^j^^ corresponds to the 
beginning of the shower as expected. This relation can 
be seen in Fig. [7] and is discussed in depth later in this 
section. 

In the analytical solution (Fig. IH]), for that particular 
observer seeing the region around shower maximum with 
angles close to the Cherenkov angle, it is evident that 
a given observation time t corresponds to two different 
shower coordinates Zj_ one corresponding to an early de- 
velopment of the shower observed at angle 9 < 9c and the 
other to a late one at angle 9 > 9c- When viewing parti- 
cle showers around the Cherenkov angle in the near-field, 
the radiation due to the early parts of the shower inter- 
feres with radiation due to the late parts of the shower. 
The apparent violation of causality is a relativistic effect 
due to the index of refraction of the medium being n > 1. 
Note also that for shower positions observed below the 



Cherenkov angle the derivative dt/dt' < meaning that 
time appears to run backwards. 

The depths z'^ are easily obtained by expressing the 
source time t' in terms of the observer time t by inverting 
Eq. (El: 



V 



clt 



(v2 



± C, 



V(z- wi)2 -r2(n2/32 _ 1) 
(u2 - cl) 

(29) 

where c„ = c/n. Real solutions exist if the argument 
of the square root is non-negative which is equivalent to 
i > tmin with ^min given in Eq. (^5]) . 

The features of the vector potential due to the appar- 
ent motion of a charge distribution along an axis with 
V > c/n are illustrated in Fig. [SI In the bottom panel 
we show the Greisen parametrized longitudinal profile of 
the excess charge as a function of z' or equivalently t'. 
In the center panel we show the relation between t and 
t' given in Eq. (pi)) for ice (n — 1.78) with the shower 
front traveling at the speed of light. This relation has 
a minimum at t^jj^ = '^m'm/''^ given in Eq. ((26)) which 
corresponds to an observer time imin given in Eq. (j28l) . 
Around t'^^^ the derivative dt/dt' is very small and as a 
consequence the charge density of the source that cor- 
responds to a relatively wide interval of source times 
t^jj^ ± Ai' = i'-cain ^ ^z' /v, is projected and seen by the 
observer during a small interval tmin i At where. 



At 
At' 



dt_ 

dt' 



(30) 



This causes a time compression which enhances the radi- 
ation, especially when the geometry is such that at tmin 
the observer located at (r, 0, z) sees the region around 
the shower maximum. This is illustrated in the right 
panel of Fig. [5] where the vector potential corresponding 
to the charge distribution shown in the bottom panel is 
depicted. The sharp initial peak of the vector potential 
is due to the compression of radiation into a short in- 
terval of time as seen by the observer. For late enough 
observer times, the relation between the apparent and 
the observer time approaches linearity and the relation 
is causal. The observer sees later parts of the shower, 
shown in the bottom panel of Fig. [51 at later times. This 
corresponds to the long tail in the vector potential de- 
picted in the right panel of Fig. ^ It is worth remarking 
that the relativistic effects described here arise in any 
situation where v > c/n and are not limited to the de- 
scription of Askaryan radiation (sBj . 



B. Interpretation of simulation results 

The corresponding time delay analysis of the results 
of the ZHS particle simulation shown in Fig. [5] is dis- 
played in Fig. [3 using the same profile and source po- 
sitions as depicted in Fig. [SJ The observer located at 
{x,y,z) = (10,0,20) m in Figs. [5] and [3 sees a sharp and 
strong spike in the vector potential (and electric field) 
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that is not matched by that seen by observers located at 
other positions in Fig. [5] In the middle panel of Fig. [7] 
we show the corresponding observer time t vs. the source 
position z' relation. For the observer at (10, 0, 20) m one 
can clearly see a region with a small derivative dt/dt' 
responsible for the compression effect which leads to an 
enhancement in the vector potential. As shown with the 
aid of Eq. (|27p the time at which an observer sees the 
shower first corresponds to observation at the Cherenkov 
angle. This can be also seen in the top panel of Fig. [7] 
where we have plotted the angle between the position z' 
along the shower axis and the location of the observer. 
An observer at (10, 0, 10) m also sees a fraction of the 
shower with angles around the Cherenkov angle. The 
Cherenkov pulse is not as pronounced because the net 
charge in that region is significantly smaller than the 
charge at shower maximum, and does not last as long 
as what the observer at (10, 0, 20) m sees. In this view 
Cherenkov radiation is a geometrical phenomenon due to 
the minimum in the relation between the observer time 
and the source positior0. 

The apparent causality violations are manifested in 
the shape of the vector potential as viewed by differ- 
ent observers. An observer at an angle smaller than the 
Cherenkov angle will see the evolution of the shower with 
an inverted causal order at all times. This is the case of 
the vector potential labeled {x,y,z) — (10,0,40) m in 
Fig. [S] where the observer sees a ^ 25 m long 100 PeV 
energy shower (shown in the inset of Fig. [5] and in the 
bottom of Fig. [7]) with the radiation corresponding to 
larger depths at later times. The longitudinal charge 
excess distribution shown in Fig. [7] has a primary peak 
followed by a smaller secondary peak. The vector po- 
tential traces this feature but in the reversed time se- 
quence. In the case of an observation at angles larger 
than 9c, the shower is seen in the normal causal order 
at all times. This case is illustrated by vector potential 
labeled (x, y, z) = (10, 0, 10) m in Fig. [5] where the radi- 
ation corresponding to the charge excess distribution at 
larger depths is observed at later times. The peaks of 
the vector potential match the order of the peaks of the 
charge excess distribution in the bottom of Fig. [T) It is 
also worth noting that even though single antenna mea- 
surements in the near-field make it difficult to reconstruct 
the longitudinal profile of the charge excess, observations 
from multiple stations located at tens of meters from the 
shower axis, do provide the necessary information. 



^ One cannot forget that the geometrical effect manifested as an 
index of refraction n > 1 is in fact due to an interaction of the 
excess charge with the atoms in the medium and which is, in 
general, frequency dependent. The discussion in this section is 
only relevant to frequency bands where the index of refraction is 
reasonably approximated by a constant. 




Source Position z (m) 

FIG. 7. Apparent motion analysis of the charge distribution 
for the vector potentials shown in Fig. [5] Bottom panel: Lon- 
gitudinal profile of the excess negative charge in the 100 PeV 
electron-induced shower. Middle panel: Observer time t vs 
source position z' for the same observers as in Fig (5] Top 
panel: Angle with which each observer sees the different posi- 
tions z' in the shower axis. The Cherenkov angle (6c ~ 55.8° 
in ice) is indicated with a long-dashed horizontal line. See 
text for a description of the mapping of these relations to the 
features of the vector potentials shown in Fig. \5\ 



IV. SUMMARY AND OUTLOOK 

We have derived a highly detailed and computationally 
efficient approach for the calculation of Askaryan pulses. 
The electrodynamic calculations leading to the relation 
between the pulse features and the shower characteristics 
can be intuitively understood via the apparent motion of 
charges. Viewed from this perspective, one can easily re- 
trace the time-domain behavior of the pulse to the shape 
of the electromagnetic shower and the observer versus 
source time relation as shown in Fig. El 

There are many interesting features of Askaryan radi- 
ation due to the fact that the speed of the shower front 
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exceeds the speed of light. In this treatment, the ra- 
diation due to the charge excess of an electromagnetic 
shower is understood as a "dense" or compressed map- 
ping of the charge excess profile to the vector potential 
via the observer vs. source time relation. The mapping 
is densest at the minimum of the observer vs. source 
time relation which corresponds to observations at the 
Cherenkov angle. In addition, for observation at an- 
gles smaller than the Cherenkov angle, time appears to 
run backwards. This manifests itself in the time-reversed 
mapping of the longitudinal profile of the charge excess 
to the time-domain vector potential. 

The primary motivation for developing this calculation 
in the time-domain is to understand the temporal be- 
havior of the Askaryan pulse. In the frequency domain, 
this is equivalent to understanding the phase versus fre- 
quency relation. Although it is possible to do this in a 
completely frequency-domain approach, the time-domain 
relations can be intuitively understood and are easier to 
compute. 

The computational algorithm presented here is sum- 
marized as follows: 

1. Compute the vector potential of the Askaryan pulse 
at the Cherenkov angle A{6c, t) and use it together 
with total charged track-length LQtot = J dz'Q{z') 
to extract the functional form of the form factor 
Fp. (This has been done in this article for elec- 
tromagnetic showers in ice and in that case it is 
possible to use directly the parameterization pro- 
vided in Eq. (fT6)) . For other situations it needs to 
be re-evaluated with a detailed Monte Carlo simu- 
lation.) 

2. Obtain the charge excess longitudinal profile of an 
electromagnetic shower Q{z'). This can be pro- 
vided as either the output of a particle shower sim- 
ulation or using a parameterization. 

3. Convolve Fp with Q{z') according to Eq. p7|) in 
the far-field or Eq. (|22l) in the near-field to obtain 
the time-domain vector potential. 

4. Electric fields are obtained from a trivial numerical 
derivative of the vector potential with respect to 
time: E = -dA/dt. 

The formalism developed here can also be applied to 
the reconstruction of longitudinal shower profiles. In 
the far-field this can only be done if the pulse was de- 
tected away from the Cherenkov angle, otherwise the 
pulse shape is approximately the same for any given lon- 
gitudinal charge excess profile. Away from the Cherenkov 
angle the pulse traces the shape of the longitudinal profile 
convolved with the lateral profile response. For very ex- 
tended longitudinal profiles, such as those resulting from 



UHE showers affected by the LPM effect, the tracing of 
the profile can be seen for angles as small as 0.3° away 
from the Cherenkov angle (see Fig. [3]). The overall qual- 
ity of the reconstruction can be assessed with simulations 
that are specific to the experiment in question. The re- 
construction of longitudinal profiles has interesting ex- 
perimental applications such as the identification of the 
primary particle or v flavor inducing the shower. This is 
particularly relevant for the electromagnetic component 
of a z/g-induced shower with its multiply-peaked structure 
due to the LPM effect. 

In the near-field the reconstruction of shower longitu- 
dinal profiles is also possible. If the radiation is detected 
by a single station the reconstruction is complicated by 
the fact that the portions of the shower above and be- 
low the Cherenkov angle may interfere with each other 
depending on the position of the antenna. However, if 
multiple antennas observe the radiation due to a single 
shower it is possible to obtain a highly constrained re- 
construction of the longitudinal profile. The example de- 
picted in Fig. [S] shows a case where this could be done for 
antennas in ice spaced tens of meters apart. The formal- 
ism provided in this paper will allow the experimentalist 
to simulate these measurements and find the optimal an- 
tenna placement and assess the quality of reconstruction. 
This would be of particular interest to an array such as 
the future planned ARA 14] and ARIANNA [ll] experi- 
ments where the shower could potentially be observed by 
multiple stations in Antarctic ice. 

In a future publication we plan to produce time- 
domain parameterizations of the vector potential at the 
Cherenkov angle for electromagnetic showers in various 
media such as salt and the lunar regolith. In addition the 
time-domain parameterizations of hadronic showers will 
be included which will allow the experimentalist to pro- 
duce full simulations of neutrino interactions with flavor 
dependent parameters. This will be useful for experiment 
simulations and candidate event reconstructions. 
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